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Abstract: 

Macroscopic quantum coherence oscillations in mesoscopic antiferromagnets 
appear when the anisotropy potential creates a barrier between the antifer- 
romagnetic states with opposite orientations of the Neel vector. This phe- 
nomenon is studied for the physical situation of the nuclear spin system of 
eight Xe àtoms arranged on a magnètic surface along a chain. The oscillation 
period is calculated as a function of the chain constant. The environmen- 
tal decoherence effects at finite temperature are accounted assuming a dipole 
coupling between the spin chain and the fluctuating magnètic field of the sur- 
face. The numerical calculations indicate that the oscillations are damped by 
a rate ~ (N — l)/r, where iV is the number of spins and r is the relaxation 
time of a single spin. 
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I. Introduction 



The observation of macroscopic quantum coherence (MQC) phenomena 
in complex many-particle systems represents a subject of wide interest, rang- 
ing from the conceptual foundations of quantum mechanics 0, to the physics 
of the microelectronic devices. During the last years a particular attention 
was given to the macroscopic quantum tunneling and quantum coherence 
oscillations [0]. In these phenomena the quantum dynamics of a prepared 
non-stationary wave-packet is directly reflected by the non-classical behavior 
of a macroscopic observable. However, it is dimcult to identify physical situ- 
ations where such coexistence of the classical and quantum aspects could be 
observed. 

Quantum coherence oscillations may occur in the localization probabil- 
ity of relatively complex systems as the individual Xe àtoms trapped in the 
surface-tip junction of the scanning tunneling microscope [§, |||, g |(J. Partic- 
ularly suitable candidates to observe MQC phenomena are also the magnètic 
systems J7|. In the anisotropic antiferromagnets the Neel vector may change 
the orientation by quantum tunneling || or quantum coherence oscillations 
§]. The observation of these oscillations still represents a challenging prob- 
em, but important results have been obtained from the measurements of the 
ac magnètic susceptibility in antiferromagnetic (AF) ferritin fpL0| 3 p . 

The study of the non-elementary excitations in ÀF systems is also essen- 
tial for understanding the high-T c superconductivity [12], |T3| . Therefore, a 



detailed study of the MQC phenomena in antiferromagnets appears highly 
interest ing. 

The purpose of this work is to investigate the occurence of MQC oscilla- 
tions of the Neel vector in an anisotropic AF chain C = {Ij, i — 1, 8} of eight 
spins 1/2. The Hamiltonian correspond to the system of nuclear spins for 
eight Xe àtoms placed on a magnètic surface and coupled by the magnètic 
dipole interaction. Such chain structures could be constructed, for instance, 
using the scanning tunneling microscope in the "atòmic switch" operation 
mode, proved during the last years to be an efficient instrument to displace 



in a controlled way the Xe àtoms on a Ni surface |L4 

The model Hamiltonian and the spectrum of the spin chain are presented 
in Sect. II. It is found that the first two levels are quasi-degenerate, and 
dominated by AF states with opposite orientations of the Neel vector. There- 
fore, non-stationary AF wave packets which are interchanged periodically by 
MQC oscillations can be constructed. Their oscillation period is calculated 
as a function of the chain constant. 

The relaxation effects at finite temperature due to the magnètic coupling 
between the spin chain and the surface are discussed in Sect. III. The cou- 
pling between the electrònic and the nuclear magnètic moments is significant, 
and in the opposite situation, when the nuclear spins play the role of the envi- 
ronment, the MQC oscillations in small magnètic partides can be suppressed 
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15]. Here the environment is represented by the surface, and the coupling is 
renected by the finite relaxation time r of the nuclear spins. The effects of 
the spin relaxation on the nuclear MQC oscillations will be accounted by the 
coupling between the chain and the fluctuating magnètic field of the mag- 
nètic surface. It is found that when r is in the range of seconds the MQC 
oscillations are damped, but not completely suppressed. The conclusions are 
summarized in Sect. IV. 



II. Quantum coherence oscillations in the AF nuclear spin chain 

The chain considered in the present calculations consists of N = 8 spins 
1/2 equally spaced by d, and interacting by the model Hamiltonian 

N-l 

Í + 1,X Cyíi^yí 

i+l,y + C z l itZ ï i+líZ } . (1) 

i=l 

Here J > is the coupling strength and c x , c y and c z are anisotropy coef- 
ficients. Classically, the AF ordering is described by using two sub-chains, 

C and C e , containing the odd, {Ij, i = 1, 3, 5, 2N — 1} and the even spins 
{li, i = 2, 4, 6, 2A/"}, respectively. The corresponding magnetization vec- 
tors are M Q = 7^Z)i=oddIj an d M e = 7^ Si= et >en where 7 denotes the 
gyromagnetic factor. The equilibrium configuration of this system is antifer- 

romagnetic when M Q and M e have the same magnitude Mq, but a relative 
antiparallel orientation. The energy in this case is a function of the Neel 

vector ïï = (M — M e )/2M , and is expressed by 

E A (n) = -J(N- l){c x n 2 x + c y n 2 y + c z n 2 z ) / 4 . (2) 

An anisotropic system with c Xty < c z has two degenerate minima, E^ in = 
—J(N — l)c a /4 attained when n has the two possible orientations along the 
Z-axis, n = ±é*2. These minima are separated by a two-dimensional potential 
barrier with the maximum E A ax ((p) = —J(N — l)(c x cos 2 4> + c y sin 2 0)/4, in 
the X-Y plane, where denotes the angle between the Neel vector and the 
X-axis. 

A physical situation of interest described by the anisotropic Hamiltonian 

of Eq. (1) appears when C = {li, i = 1,8} consists of the nuclear spins for 
a chain of Xe àtoms placed on a magnètic surface. Such chains can be con- 
structed, for instance, using techniques of atòmic manipulations. The inter- 
action between the Xe nuclear spins is due to the magnètic dipole forces, and 
the effect of the surface can be taken into account using the simple method 
of images. In a coordinate system with the Z-axis normal to the surface and 
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the Y-axis along the chain, each magnètic moment m = (m x ,m y ,mz) lying 

above the surface at the distance z = tq = 2.17 À, equal to the radius of 
the Xe atom, has an image m! = (—m x , —m y ,m z ) at z — —tq. The nearest 
neighbors of m interact also with mf, and the Hamiltonian of the whole sys- 
tem is expressed by Eq. (1), where J = (J>oh 2/ y 2 /(Attcí 3 ), fio = 4irl0~ 7 N/A 2 is 
the vacuum permeability, 7 = —I-òA/in/H is the gyromagnetic factor of the 
129 Xe isotope and d is the chain constant. The anisotropy coeficients are 

c x = 1 — 2 sin 3 a 
c y = -2 - 2sin 3 a(l - 3sin 2 a) (3) 
c z = 1 + 2 sin 3 a(l — 3 cos 2 a) 

with tana = d/2r . 

The spectrum of H was calculated for 6.8 À< d < 9.4 À by solving 
numerically the eigenvalue equation 

H \if> n >= E n \ÍJ n > . (4) 

The basis was defined by the common eigenstates for the Z-components 
of all the spin operators Ij )2 , i=l,8. These basis states are denoted by 
\k >= \rrii, 7712, m*, m'l, m^, m|, m*, m% >, k — 1,2 8 , with m\ = ±1/2. 
The eigenstates of Hq have the general form \ip n >= J2kUk\^ >• ^ i s found 
that the ground (n = 0), and the first excited state (n = 1), have the largest 
overlap with only two basis states, denoted | f> and | |>, which are antifer- 
romagnetic in the classical sense. Explicitly these basis states are 

.1 1 1 1 1 1 1 1 , 

T>-| 2' 2'2' 2'2' 2'2' 2 > ' ^ } 



and 



1 1 _1 1 1 1 1 1 

i>- I - 2' 2' _ 2' 2' ~2' 2' ~2' 2 > ' 



and are eigenstates of Z-component of the Neel operator N = X^=i,4(l2i-i — 

I 2 j)/4 with eigenvalues +1 and —1. 

The results will be presented in detail for the particular case of d = 7 
À, when J = 0.17^/ms and the anisotropy coeficients given by Eq. (3) are 
c x = —0.22, c y = —0.58 and c z = 1.2. The overlap coeficients between 
the eigenstates \ipo > and > and the AF states of Eq. (5) and (6) are 
<| |^ >=<| >= 0.57, - <T |^i >=<| >= 0.63. Thus, the two 
AF states exhaust more than 65 % of the eigenstates norm. The eigenvalues 
E = — 0.406/i/ms and E 1 = — 0.395/i/ms are separated by A = E\ — E — 
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O.OlUt/ms which is sensibly smaller than E2 — E\ = 0.053/ï/ms. The small 
value of A shows that the system has a quasi-degenerate ground state. This 
appears as a "tunneling doublet" determined by the two-dimensional poten- 
tial barrier E max ((f)) = (0.066 cos 2 4> + 0.17sin 2 <f))h/ms separating the AF 
energy minima of E min = — 0.353/ï/ms. 

The tunneling behavior is shown clearly by the evolution of the non- 
stationary wave-packets prepared at the ground-state energy with a well- 
defined AF configuration. States with these properties are represented by 
the linear combinations 

>= —^(1-00 > +|^i >) , |Vr >= ^(IV'o > —1-01 >) • (7) 

of the eigenstates \ipo > and j^i >. These wave-packets can be distin- 
guished macroscopically by the expectation value of N z , and during the 
time-evolution are interchanged periodically, by MQC oscillations. Thus, 
if the system is prepared at t=0 in the state >, then at the moment t it 
will be found in the state |V>t > with the probability V^{t) = sm 2 (7it/2T max ), 
where T max = hn/A. The half-period of oscillation T max is represented as 
a function of the chain constant d in Fig. 1 (A). For d = 7 A, one obtains 
T max = 0.3 s, and the calculation of the expectation vàlues 

< N k > (t) =< ^ í \é Hotlh N k e- iHotlh \^ l > (8) 

shows that in time < N x > (t) = 0, < N y > (t) = 0, while the Z- 
component has MQC oscillations < N z > (i) =< ipo\N z \ipi > cos(7it/T max ), 

with < Vp|AyV>i >= ~ - 83 ( Fi S- 1 ( B ))- 

It is important to emphasize the extreme sensitivity of the MQC res- 

onance with respect to the preparation of the initial state. The expected 
value of the energy in the classical antiferromagnetic states, <| \H \ |> and 
<t \H \ t>, is E min = -0.353^/ms, higher than E x = -0.39^/ms by w 4A. 
This energy is far outside the interval [Eq,Ei], and therefore between the 
classical antiferromagnetic states | |> and | f> there are no MQC oscilla- 
tions. 



III. MQC oscillations at finite temperature 



The evolution of the Neel vector at the MQC resonance, presented in the 
previous section, was obtained neglecting the residual interactions between 
the nuclear spins and the environment. However, when these interactions 
are considered the MQC oscillations could be damped |T^] or completely 
suppressed |T^|, and the Neel vector aquires a fixed orientation. At finite 

temperatures the nuclear spin Ij of each Xe atom lying on the magnètic 
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surface is subject also to an external, time-dependent magnètic field Bf(t), 
created by the phonon modulation of the crystalline elèctric field and the lat- 
tice spin waves |R| . Due to this field, appears in the Hamiltonian a residual 
interaction term 

flr(t)=-7»EÏ^(í) • (9) 

i=l,8 

At low temperatures the typical wave length of the surface phonons and 
magnons is ~ 400 À, much larger than the length of the chain, L = 49 

A. Therefore will be considered to be the same for all spins, Bf(t) = B e (t). 
With these approximations the residual interaction term determined by the 

environmental magnètic field is H r (t) = — r yKB e (t) ■ I, where I = J2í=i s !«• 

The "atòmic switch" experiments [P[ have been performed at the envi- 
ronmental temperature T = 4K, when the thermal energy /c^T = 0.34 meV is 
very high compared both to the tunnel splitting A = 6.83 ■ 10~ 12 meV and to 
the maximum barrier height Vb = E^ ax (7i/2) = 0.11 • 10 -9 meV. Therefore, 
the thermal environment can be considered as classical, and the field compo- 
nents B e {t), [i — x, y, z, will be treated as a white noise with zero mean. The 
normalization of this fluctuating field is ensured by the fluctuation-dissipation 
theorem (FDT) « B*(t)B*,tf) »= 5^8{t - t')/(-y 2 T), where « ... » 
denotes the average over the statistical ensemble describing the environment. 

For a single spin of the chain, H r (t) induces transitions between the states 
|l/2 > and | - 1/2 > with a rate |j À = (| < -l/2|I x |l/2 > | 2 + | < 
—l/2\ï y \l/2 > | 2 )/r = l/2r, and the relaxation rate of the population differ- 
ence n\ii — n-1/2 is 2À = 1/r. Therefore, the parameter r has the meaning 
of spin-surface relaxation time, increasing as 1/T when T — > 0. 

The relaxation of the MQC oscillations cannot be treated by using a sim- 
ilar two-level approximation, because the matrix elements of H r within the 
subspace generated by | Vo > an d l^i > are 0. The operators ï x and l y con- 
tained in H r act on the initial state > by flipping the individual spins, 
and in time the quantum state aquires components over the whole spectrum. 
This process is described in principle by a transport equation for the den- 
sity matrix 0, |2Ü| but due to the relatively large number of states (= 2 8 ), 
such numerical calculations are not feasible. Moreover, at high temperatures, 
when the initial state of the system is a pure state and there are only few ob- 
servable of interest, as in the present case, the computational effort required 
to find the evolution of the whole density matrix is not justified. Instead, an 
equivalent description pi which can be applied efficiently is provided by a 



statistical ensemble of N t Brownian trajectòries \il> r (t) >, r = 1, N t , obtained 
by integrating the Schrodinger equation 

ihd t W{t)>=[H ü - 1 %B e {t)-ï]\i) r {t)> . (10) 
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The numerical integration was performed using the procedure presented in 
[H], as a classical system of Hamilton equations for the real and immaginary 
parts of the amplitudes Vk(t) =< k\if)(t) >, k = 1,2 8 . Using the notation 
ttfc(i) = Re(yk(t)) and v k (t) = Im(y k (t)), Eq. (10) takes the form 

2hu k = — 2hv k = , (11) 

avk ou k 

where 

W(í) = E i u ^k' + v k v y )Re{< k\H + H r {t)\k' >)- 

k,k'=l 

{u k v k > - v k u k >)Im(< k\H + H r .(t)\k' >) . (12) 

The dominant AF states are annihilated by I 2 , and therefore the contribution 
of the term B e z (t)\ z from H r is very small, and it was neglected. The remain- 
ing X and Y components of the nuctuating field at the moment t n = ndt, 

normalized according to the FDT, have the form B e {t n ) = R n \j ! \j '(j 2 rdt) 

where {Rn, n = 1,2,3,....} is a sequence of random numbers with mean 
and variance 1. 

The time-evolution of the ensemble average for any observable O is de- 
fined by 

i m 

«0»(t) = -J2<^ r (t)\0\^(t)> . (13) 

iV í r=l 

The observables of interest here are the Neel vector and the energy, and their 
averages were calculated using Eq. (13) with O denoting the components of 
the Neel operator N k , and H , respectively. The results obtained for d = 7 
A, r = 2.5 s and N t = 20 are presented in Fig. (2). The components 
<< N x >> and << N y » of Fig. 2 (A) and (B) have thermal fluctuations 
around 0, while << N z » shown in Fig. 2 (C) has damped oscillations. 
These oscillations can be well approximated by the analytical expression 

« N z » (t) = e~ At < N z > (t) , (14) 

where the damping constant obtained by fit is A = 3 s -1 . Similar calculations 
with t in the range of seconds, as expected at low temperatures, indicate that 
A = 7.5/r. 

The ratio (<< H » (t) — E )/A between the average excitation energy 
and the doublet splitting À is a suitable measure of the heating effect pro- 
duced by the environment, and is represented in Fig. 2 (D). The ensemble 
average of the energy appearing here is accurately reproduced by the formula 
<< H » (i) =<< H » (0) + A(wit - w 2 t 2 ), with w\ = 45.9 s" 1 and 
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W2 = 19.3 s~ 2 . 



IV. Conclusions 

In this work it was shown that the chain of nuclear spins for eight àtoms 
of the 129 Xe isotope placed on a magnètic surface has a low-energy mode of 
collective excitation represented by antiferromagnetic MQC oscillations. The 
numerical calculations indicate that when the chain constant d is larger than 
6.8 A the flrst excited state ipi of the chain has a very low excitation energy. 
This state is strongly coupled to the ground state tpo by the Z-component 
of the Neel operator, such that | < ipo\N z \ip 1 > \ = 0.83. Both states, ipo 
and ipi, are dominated by two AF spin configurations having opposite ori- 
entations of the Neel vector. Therefore, resonant wave packets containing 
only one of these AF configurations can be constructed by the symmetric 
and antisymmetric linear combinations of ipo and ipi. These wave-packets 
are non-stationary, and are interchanged periodically in time. Although the 
number of spins considered here is relatively small, the resonant wave packets 
are distinguished by the orientation of the Neel vector, and their oscillation 
represent a MQC phenomenon. When d = 7 A the half-period of oscillation 
is T max = 0.3 s, and increases exponentially with d (Fig. 1 (A)). 

The MQC oscillations are known to be very sensitive to the decoherence 
effects produced by the coupling to the thermal environment. If in practice 
these effects are strong, then the oscillations may be completely suppressed 
and cannot be observed. In the present case decoherence appears due to 
the residual dipole interaction between the nuclear spins and the fluctuating 
magnètic field produced by the surface. This coupling was described by one 
parameter, chosen to be the relaxation time r of a single nuclear spin. 

The relaxation rate of the MQC oscillations can be easily calculated in the 
physical situations which can be treated within a two-level approximation. 
However, for eight spins the number of states coupled by the residual dipole 
interaction term is large, and this approximation cannot be applied. There- 
fore, the fluctuating residual interaction term was included in the Hamilto- 
nian, and the diffusion of the initial wave-packet \ipi > was described by 
an ensemble of Brownian solutions obtained by integrating the Schròdinger 
equation. The results indicate that when r is in the range of seconds, the 
MQC oscillations are not suppressed, but are damped by a rate ~ (iV — l)/r, 
where N is the number of spins. 
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Figure Captions 

Fig. 1. The half-period of the MQC oscillations T max as a function of the 
chain constant d (A) and the expected value < N z > as a function of time 
when d = 7 À (B) . 

Fig. 2. Ensemble averages of the Neel vector components << N x >> 
(A), << N y » (B), << N z » (C) and (« H » -E )/A (D) as 

functions of time when d = 7 À and r = 2.5 s. 
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